Hysteretic transition from laminar to 
vortex shedding flow in soap films 

V.K. Horvath*, J.R. Cressman, W.I. Goldburg, and X.L. Wu 

Department of Physics and Astr., Univ. of Pittsburgh, Pittsburgh, PA 15260, USA 
"Contact e-mail: vhorvath@pitt.edu 

1 Introduction 

Studies of flow behind a single cylinder have greatly contributed to our under- 
standing of the development of complex flows like the Karman vortex street jlj 
or turbulent flows. At the first instability in such a system the laminar flow (LF) 
becomes time dependent and vortices appear. It has been shown that this tran- 
sition to the vortex shedding phase (VS) is well described by the supercritical 
Hopf bifurcation Q]. 

2 Hysteretic Primary Instability in Soap Films 

Soap films have been shown to be a particularly useful for the study of two dimen- 
sional flows H). Here we present measurements demonstrating that the LF— >VS 
transition can be hysteretic in a (quasi-two dimensional) soap film penetrated 
by a glass cylinder [|J. Our experimental setup consists of a rapidly flowing soap 
film formed between two vertically positioned thin nylon lines (sec Fig. 1). The 
film is fed continuously from the top through a valve. To investigate the flow 
at different Reynolds numbers Re = Vd/v we changed the mean velocity V by 
cither opening and closing the valve at different rates or by changing the separa- 
tion of the nylon lines, using computer controlled stepping motors (here v is the 
kinematic viscosity). At the center of the parallel segment of the set-up, a glass 
rod penetrates the film in the z direction. The velocity measurements were made 
using a dual head laser Doppler vclocimctcr (LDV). The fluctuating horizontal 
component V x (t) and V were measured simultaneously below and above the rod 
respectively. 

On increasing the flow rate a transition appears from the LF state to the 
VS state. Below the transition, the flow is characterized by two counter-rotating 
vortices underneath the rod (Fig. 2a). In the VS state these recirculating vortices 
peel off the rod and flow downstream. The continuously generated counter- 
rotating vortices form the Karman vortex street (Fig. 2b). All experiments start 
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Figure 2: Interference pictures of the 
Figure 1: Schematic of the experimen- flowing soap film below the rod in (a) 
tal setup. The diameter d of the rod is the laminar state and (b) the vortex 
1mm. The 0.2mm thick nylon lines are shedding state. 

w?l^slow^flOTcasmg ^^m^hc LF regime. The critical velocity, where VS 
commences, is called V^ p . After some waiting time in the VS state, V is slowly 
decreased. At V^ own the system undergoes a reverse transition, i.e. from VS 
into LF. This cycle was repeated several times in each run with V x {t) and V 
being recorded simultaneously. 

In order to determine the transition velocities, we have calculated the velocity 
probability distribution function P(V X ) from V x (t) by a standard binning proce- 
dure. In Fig. 3 the magnitude of P(V x ,t) is mapped into gray scale values and 
is shown as a function of V x and t. It can be seen that sharp changes in P(V X , t) 
precisely indicate the transition times t* d and f* and therefore V^ own = V{t* d ) 
and V ( up = V(t*) too. It is apparent in Fig. 3 that V c up is not equal to V c down , 
that is, there is no unique critical velocity for the transition. The Reynolds 
number is roughly 50 in the hysteretic gap. A large uncertainty in this value 
comes from the poorly determined two-dimensional soap film viscosity v , which 
depends on the film thickness. 

It is worth pointing out that if we keep the mean flow rate constant in the in- 
terval [Vf own ,V ( ! J ' p ], then the laminar state can persist for an indefinite length of 
time. However applying sufficiently large acoustic or pulse-like mechanical per- 
turbations, the system could be driven into the VS state. The system remains 
in this state even if the perturbation is turned off. Applying similar pertur- 
bations, VS-^LF transition was never observed at constant V. As a result we 
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Figure 3: Here the function V(t) is 
shown together with the P(v x ,t) map 
at 1^/^1=0.004 m/s 2 . The gray scale 
corresponds to the probability density 
of v x . The solid and dashed lines are 
intended to guide the eye for the tran- 
sitions VS^LF and LF^VS. 
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Figure 4: This figure shows our experi- 
mental data together with different fits 
(see text). The circles and squares rep- 
resent two different experimental run. 
Closed and opened symbols are used for 
data taken at increasing and decreasing 
V respectively. 



can conclude that the observed static hysteresis is not a result of some delay in 
the response of the system to the bifurcation. Auxiliary experiments were also 
performed to exclude wetting properties and the effect of the air boundary near 
the film as possible sources of hysteresis. 

In the absence of any available theory, we consider a fifth order Landau 
equation^ dv(t)/dt = av + 2bv 2 v — cv 4 v, for the complex velocity v~v x + iv y . 
Here all variables denoted by tilde are complex and — Vi — < Vi >t (i—x or y). 
The complex velocity can also be written as v=ve v ? ', where v is the amplitude. 
The real part of this equation describes the evolution of v. Its non-trivial positive 
stationary solution describes the VS state: 

VVS 2 = - (l + ^l+a r C r /bl, (1) 
C r 

where the subscript r is used to designate the real part of complex numbers. 

Although the application of this model for data fitting purposes is not obvious 
(the relationship between the parameters a r ,b r , c r and the control parameter of 
the experiment is unknown), typically a r is considered to be proportional to 
(V — V™ p ). As one can see in Fig. 4, this generic model of the hysteresis is 
in qualitative agreement with our observations. However none of our 3 best 
fits are very satisfactory. It is interesting to note that a simple three-parameter 
phcnomcnological equation Pq (V — V c ) + 1P\v — P2V 2 = provides a surprisingly 
good fit to our data in the VS state (see inset of Fig. 4.). Using least square 
fitting to the solution of this equation provides the following result: v x = 0.04 ± 

yV/39.06 — 0.0113, where the velocities are in units of m/s. 
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Auxiliary experiments suggest that the hysteresis may be connected with the 
fact that the film may become very thin in the recirculating region. Recognizing 
that the viscosity of a soap film depends on its thickness, the observed hysteresis 
effect can be qualitatively understood; the Reynolds number can now take dif- 
ferent values at the same value of V, when the parameter is lowered than when 
it is raised. 

3 Conclusions 

• We have observed unexpected hysteresis at the onset of the Karman vortex 
street in a quasi two-dimensional soap film. 

• The fifth order amplitude equation is in qualitative agreement with our 
observations, but a simple phenomenological equation provides a better 
numerical fit to the experimental data in the vortex shedding regime. 

• A phenomenological picture is suggested to explain the origin of the hys- 
teresis. 
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